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Abstract.  Cyclic  signal  processing  refers  to  situations  where  all  the  time  indices  are  interpreted  modulo  some 
integer  L.  In  such  cases  the  frequency  domain  is  defined  as  a  uniform  discrete  grid  (as  in  L-point  DFT). 
This  offers  more  freedom  in  theoretical  as  well  as  design  aspects.  While  circular  convolution  has  been  the 
centerpiece  of  many  algorithms  in  signal  processing  for  decades,  such  freedom,  especially  from  the  viewpoint 
of  linear  system  theory,  has  not  been  studied  in  the  past.  In  this  paper  we  introduce  the  fundamentals 
of  cyclic  multirate  systems  and  filter  banks,  presenting  several  important  differences  between  the  cyclic 
and  noncyclic  cases.  Cyclic  systems  with  allpass  and  paraunitary  properties  are  studied.  The  paraunitary 
interpolation  problem  is  introduced  and  it  is  shown  that  the  interpolation  does  not  always  succeed.  State 
space  descriptions  of  cyclic  LTI  systems  are  introduced,  and  the  notions  of  reachability  and  observability 
of  state  equations  revisited.  It  is  shown  that  unlike  in  traditional  linear  systems,  these  two  notions  are  not 
related  to  the  system  minimality  in  a  simple  way.  Throughout  the  paper,  a  number  of  open  problems  are 
pointed  out  from  the  perspective  of  the  signal  processor  as  well  as  the  system  theorist. 


1.  INTRODUCTION 

Consider  two  sequences  x{n)  and  h{n)  defined  for  0  <  n  <  L  —  1,  and  let  j/(n)  denote  their  circular  or  cyclic 
convolution  [17].  That  is,  y{n)  =  X]m=t)a;(m)/i(n  -  m)  with  all  arguments  interpreted  modulo  L.  We  can 
represent  this  operation  by  the  block  diagram  of  Fig.  1  where  x(u)  and  y{n)  are  the  input  and  output, 
respectively,  of  a  linear  system.  W^ith  all  time-arguments  interpreted  modulo-L,  this  is  also  a  time-invariant 
system  (i.e.,  a  circular-shift  invariant  system).  We  say  that  {h(n)}  is  a  cyclic  LTI  system.  For  the  purpose  of 
interpretation  one  can  also  regard  x(n)  to  be  a  periodic-L  input,  for  which  the  LTI  system  h(ji)  yields  the 
periodic-L  output  y{n). 

Circular  or  cyclic  convolution  has  been  at  the  center  of  digital  signal  processing  from  its  early  days  [8]. 
Indeed,  fast  algorithms  for  ordinary  convolution  routinely  convert  the  problem  into  a  cyclic  convolution  and 
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then  use  the  FFT.  The  notion  of  polynomial  transforms  introduced  by  Nussbaumer  [16]  is  another  example 
of  the  prevalence  of  these  ideas  in  the  early  days  of  digital  signal  processing.  Most  of  the  classical  applica¬ 
tions  involving  cyclic  convolutions  are  based  on  the  important  result  that  cyclic  convolution  corresponds  to 
multiplication  of  the  DFT  coefficients,  that  is,  Y{k)  =  H{k)X{k)  for  0  <  fc  <  L  -  1.  Farther  properties  of 
cyclic  LTI  systems  from  the  viewpoint  of  multirate  signal  processing  (more  generally  linear  system  theory) 
have  not  been  studied  in  the  past.  In  this  paper  we  wUl  consider  a  number  of  such  properties.  The  emphasis 
wUl  be  on  well-known  topics  [1],  [12],  [23],  [27]  such  as  multirate  systems,  filter  banks,  paraunitary  matrices, 
and  state  space  representations,  but  this  time  in  a  cyclic  setting. 

1.1.  Motivation  and  Scope 

The  frequency  response  of  a  cyclic  LTI  system  is  the  L-point  DFT  of  the  impulse  response. 

L-l 

H{k)  =  h{n)Wt,  0<k<L-l,  (1) 

n=0 

where  Wl  =  This  is  equivalent  to  sampling  the  conventional  frequency  response 

at  L  discrete  values  of  the  frequency  Wk  =  (the  DFT-frequencies,  Fig.  2).  The  basic  building  blocks 

in  the  implementation  of  this  system  are  multipliers,  adders,  and  cyclic  delays  Wj^.  The  cyclic  delay,  indicated 
in  Fig.  3,  has  input-output  relation  y{n)  =  x{n  -  1)  where  the  time-arguments  are  interpreted  modulo-L. 
Fig.  4  shows  the  direct-form  structure  for  H{k),  using  these  building  blocks. 

In  cyclic  signal  processing,  the  frequency  domain  is  defined  as  a  set  of  L/  discrete  frequencies  rather  than 
the  entire  range  0  <  u  <  2tt.  This  offers  more  freedom  in  theoretical  developments.  For  example,  we  will 
see  in  Sec.  4  that  the  definitions  of  allpass  filters  and  paraunitary  matrices  are  less  restricted  in  the  cyclic 
world.  Similarly  orthonormal  filter  banks  in  the  cyclic  world  are  less  restrictive.  In  subband  and  transform 
coding  problems,  if  the  power  spectrum  of  the  input  signal  is  defined  only  over  a  discrete  grid  of  frequencies 
and  the  filters  optimized  for  these  frequencies,  it  might  offer  an  increased  coding  efficiency.  Likewise,  in 
filter  design  problems  one  traditionally  imposes  constraints  (e.g.,  linearity  of  phase)  and  optimizes  the  filter 
with  respect  to  some  criterion  (e.g.,  minimax,  minimum  error  energy,  etc.).  These  constraints  are  now  only 
over  a  discrete  grid  rather  than  a  continuous  range  of  frequencies.  These  are  some  of  the  motivations  for 
considering  cyclic  signal  processing  as  a  theoretical  discipline  by  itself.  The  practical  advantages  obtained 
using  this  viewpoint  can  be  significant,  but  it  requires  further  detailed  work  to  quantify  these  advantages 
for  specific  applications.  A  unique  example  of  the  use  of  cyclic  filter  banks  (brought  to  our  attention  by  one 
of  the  reviewers)  is  the  method  of  autoregressive  spectral  estimation  in  subbands,  advanced  by  Nishikawa  et 
al.  in  1993  [14].  Our  emphasis  in  this  paper  will  only  be  on  the  theoretical  differences  between  cyclic  and 
noncyclic  systems,  especially  pertaining  to  multirate  systems,  filter  banks,  and  LTI  system  theory. 
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Related  literature.  Since  circular  convolution  is  an  integral  part  of  DSP,  we  can  regard  cyclic  LTI  filtering 
as  one  of  the  earliest  known  DSP  techniques  [8].  The  use  of  circular  filtering  in  image  subband  coding  (where 
images  have  finite  size  [7])  was  motivated  by  the  work  of  Smith  and  Sddins  on  symmetric  extension  techniques 
[19].  Circular  convolution  for  subband  coding  was  also  explicitly  considered  by  Kiya,  et  al  [11].  In  1993, 
Caire  et  al  [5]  introduced  wavelet  transforms  associated  with  cyclic  groups.  Cyclic  versions  of  two-channel 
filter-bank  orthonormality  and  power  complementarity  can  already  be  found  in  that  paper.  More  recently 
cylic  multirate  system  basics  were  explicitly  formulated  in  two  independent  conference  papers  [24],  [3].  While 
both  papers  start  with  a  common  basic  theme,  the  work  by  Bopardikar,  et  al  [3],  [4]  eventually  focuses  on 
the  wavelet  aspect  whereas  [24]  focuses  on  system  theoretic  aspects,  factorizations  and  so  forth.  The  very 
recent  work  in  [13]  focuses  on  interesting  details  of  the  two  channel  linear  phase  orthonormal  case.  In  this 
paper  the  emphasis  will  be  more  along  the  lines  of  [24-26]. 

1.2.  Notations  and  Abbreviations 

•  Wl  =  ^  with  the  subscript  deleted  if  obvious. 

•  The  integer  k  is  reserved  for  the  frequency  index  in  the  DFT  expressions.  Throughout  the  paper,  k  is 
therefore  analogous  to  frequency. 

•  In  the  figures  Wj^  denotes  one  unit  of  cyclic  delay,  that  is,  y{n)  =  x(n  —  1)  with  n  —  1  interpreted  modulo 
L.  This  is  analogous  to  or  in  standard  DSP  block  diagrams. 

•  The  abbreviation  siso  stands  for  single-input  single-output,  and  mimo  for  multi-input  multi-output. 

•  The  abbreviation  PU  stands  for  paraunitary. 

•  Matrices  and  vectors  are  denoted  by  bold  letters.  The  notations  A^,  and  A^  denote,  respectively, 
the  transpose,  the  conjugate,  and  the  transpose-conjugate  of  A. 

•  The  tilde  notation  is  defined  as  follows:  H(^)=Ht(l/z*). 

1.3.  Outline 

Sec.  2  introduces  the  basic  building  blocks  for  cyclic  DSP.  This  includes  filtering  structures,  cyclic  difference 
equations,  decimators,  expanders,  polyphase  representations  of  cyclic  filter  banks,  Nyquist  property,  and  so 
forth.  Sec.  3  considers  cyclic  versions  of  allpass  and  paraunitary  properties,  and  introduces  cyclic  orthonor¬ 
mal  filter  banks.  Sec.  4  studies  some  basic  differences  between  the  cyclic  and  noncyclic  cases.  For  example  we 
show  that  a  cyclic-LTI  system  can  be  paraunitary  or  allpass,  even  though  the  noncyclic  counterpart  (defined 
therein)  may  not  have  this  property.  This  shows  that  such  properties  are  less  restrictive  in  the  cyclic  case.  In 
Sec.  5  we  consider  the  paxaunitary  interpolation  problem.  We  show  again  that  noncyclic  interpolants  (more 
general  than  noncyclic  counterparts)  do  not  in  general  share  the  property  of  the  underlying  cyclic  system. 
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For  example  we  will  show  that  there  are  cyclic  paraunitary  matrices  which  do  not  have  FIR  paraunitary  in- 
terpolants,  though  HR  paraunitary  interpolants  always  exist.  In  Sec.  6  we  introduce  state  space  descriptions 
of  cyclic  LTI  systems.  We  also  revisit  the  traditional  notions  of  reachability  and  observability  in  the  context 
of  state  space  descriptions.  We  show  that  unlike  in  noncyclic  systems  these  concepts  do  not  have  a  simple 
relation  to  the  so-called  minimality  of  the  structure.  Throughout  the  paper  we  will  point  out  a  number  of 
open  problems  pertaining  to  cyclic  DSP  systems. 

2.  BASICS  OF  CYCLIC  DIGITAL  FILTERS  AND  MULTIRATE  SYSTEMS 


2.1.  Filtering  Structures  For  Cyclic  Digital  Filters 

Using  the  idea  that  Wj^  represents  a  cyclic(L)  delay  (analogous  to  ^  we  can  draw  structures  for  cyclic- 
LTI  systems  as  demonstrated  in  Fig.  4.  In  general  this  requires  L  multipliers.  By  expressing  the  frequency 
response  in  rational  form  we  can  sometimes  obtain  more  efficient  implementations.  Thus,  consider  the 


example  of  a  cyclic-L  transfer  function 


If  I'M  - 


(2) 


The  direct-form  structure  for  this  is  shown  in  Fig.  5,  and  the  input  and  output  of  this  system  are  constrained 
by  (1  -  bWl)Y{k)  =  (oo  -t-  d  W|)X(fc).  By  taking  the  inverse  DFT  of  this  equation  we  obtain  the  recursive 
difference  equation  (d.e.) 

y{n)  =  by{n  —  1)  -I-  aox(n)  -i-  aix(n  —  1)  (3) 


Since  the  time-indices  are  interpreted  modulo  L,  this  is  a  cyclic  difference  equation.  It  is  therefore  tricky  to 
establish  an  initial  condition  for  this  equation.  To  demonstrate  why  this  is  the  case,  let  L  =  3.  Repeated 
use  of  the  difference  equation  (and  using  the  facts  that  y{2)  =  j/(— 1),  y(3)  =  y(0)  and  so  forth)  yields  the 
conclusion 

(1  -  b^)y{0)  =  002^(0)  +  Q:i2;(l)  -I-  Q;2x(2) 


for  some  constants  d.  Thus  the  initial  condition  y(0)  is  not  arbitrary.  It  is  uniquely  determined  as  long  as 
b^  ^  1,  that  is,  as  long  as  6  ^  Wj  for  any  i.  This  condition  on  6  is  equivalent  to  the  obvious  requirement 
that  the  denominator  in  (2)  does  not  become  zero  for  any  k.  We  can  then  write  Y{k)  =  H{k)X{k),  and  the 
inverse  DFT  y{n)  is  fully  determined  for  all  n.  More  generally,  consider  the  cyclic(L)  transfer  function 


H{k)  = 


i:Lo<^nWt 


(4) 


With  the  implicit  assumption  that  H{k)  is  defined  for  all  k  (i.e.,  the  preceding  denominator  does  not  vanish 
for  any  k),  the  output  is  fully  determined  by  the  input.  In  particular  the  “initial”  condition  is  predetermined, 
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rather  than  arbitrary.  Fig.  6  shows  the  direct-form  implementation  of  this  system.  As  a  generalization  of 
the  difference  equation  idea,  we  will  study  in  Sec.  6  state  space  descriptions  of  cyclic  LTI  systems. 

Even  though  the  “initial”  condition  has  to  be  computed  separately,  the  use  of  a  recursive  structure 
often  results  in  reduced  computation.  For  example,  consider  the  cyclic  LTI  system  with  frequency  response 
This  can  be  implemented  as  shown  in  Fig.  4,  requiring  L  multipliers  and  L  -  1 
adders.  However,  we  can  find  a  more  efficient  recursive  implementation  by  rewriting 


1,-1 

H{k)  =  '^a^Wjr  = 

n=0 


l-a^ 


(5) 


The  expression  on  the  right  hand  side  yields  a  recursive  implementation  which  requires  only  two  multipliers 
and  one  adder. 


2.2.  Cyclic  Decimators  and  Expanders 

The  cyclic  decimator,  denoted  by  i  M  in  Fig.  7(a),  has  the  input-output  relation  y{n)  =  x{Mn).  We  assume 
throughout  the  paper  that  M  is  a  factor  of  L,  that  is,^ 


L  =  MK,  K  =  integer 


(6) 


With  x(n)  regarded  as  cyclic(L),  the  output  y(n)  is  cyclic(A’)  as  demonstrated  in  Fig.  7(b),  7(c).  Let  X{k) 
denote  the  L-point  DFT  of  x(n)  and  Yik)  the  AT-point  DFT  of  y{n),  that  is,  X{k)  =  'Z,nZl  x{n)W£'=,  0  < 
A:  <  L  -  1  and  Y{k)  =  0  <  fc  <  AT  -  1.  It  can  then  be  verified  (Appendix  A)  that 

M-l 

Y{k)  =  ^  (cyclic  decimator),  (7) 

1=0 

fox  0  <  k  <  K  -I,  Thus  the  DFT  Y{k)  is  obtained  by  superposing  M  shifted  copies  of  the  DFT  X{k)  where 
the  shifts  are  in  multiplies  of  K.  This  is  the  counterpart  of  the  traditional  aliasing  formula  for  decimators 
[23].  The  cyclic  expander,  denoted  by  t  iii  Fig.  8(a),  has  a  periodic-!^  input  x{n)  and  periodic-L  output 
y{n)  related  by 

^  ;  x{n/M),  n  =  mul.  of  M  (g) 

^  '  lO  otherwise. 

This  is  demonstrated  in  Fig.  8(b),  8(c).  The  corresponding  DFT  relation  is 

Y{k)  =  X{k)  (cyclic  expander),  (9) 

t  If  this  is  not  the  case,  then  decimation  by  M  would  retain  more  samples  than  one-out-of  M.  For 
example,  if  M  and  L  have  no  common  factors,  then  there  is  no  loss  of  samples  at  all,  and  the  decimated 
output  is  a  permuted  version  of  the  input  samples. 
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for  0  <  ifc  <  L  -  1.  (This  is  similar  to  the  periodic  extension  idea  in  DFT  theory  [8]).  Since  X{k)  are  the 
A'-point  DFT  coefficients  of  x(n),  the  L-point  DFT  Y{k)  has  the  period  K  -  L/M.  We  can  define  fractional 
decimators  in  a  manner  analogous  to  the  noncyclic  case  (for  example,  as  in  Fig.  4.1-10  of  [23]).  We  leave  it 
to  the  reader  to  figure  out  the  details. 


2.3.  Polyphase  Representation 

The  L-point  DFT  of  a  cyclic(L)  impulse  response  h{n)  can  be  expressed  as 

L-l  M-l  K-l  M-l 

H{k)  =  Y,h{n)Wf  =  ^ 

n=0  m=0  n-0  m=0 


for  0  <  fc  <  L  1,  where 


AT-l 

Em{k)='^h{Mn  +  m)W!^^,  Q<k<L-l. 

n=0 


Thus,  Em{k)  is  the  /('-point  DFT  of  the  mth  polyphase  component 


(10) 


(11) 


em{n)=h{Mn  +  m). 


(12) 


From  the  definition  of  Em{k)  we  see  that  it  is  cyclic(/C).  Eq.  (10)  is  analogous  to  the  traditional  Type 
1  polyphase  decomposition  H{z)  =  Similarly  the  Type  2  polyphase  form  is  given  by 

H{k)  =  WE’^'^Emik).  Thus  the  decimation  filter  of  Fig.  9(a)  can  be  redrawn  in  polyphase  form  as 

shown  in  Fig.  9(b).  The  form  in  Fig.  9(c)  is  obtained  by  using  noble  identities.  This  is  explained  next, 
including  the  reason  for  using  the  same  notation  Em(k)  before  and  after  the  decimator. 

Noble  identities.  Since  the  poljrphase  components  have  a  smaller  period  K ,  they  can  be  relocated  to  the 
right  of  the  decimators  (similar  to  the  use  of  noble  identities  in  traditional  noncyclic  case).  Thus  consider 
Fig.  10(a)  where  we  have  denoted  the  outputs  of  Em{k)  and  the  decimator  by  u(n)  and  t;(n)  respectively. 
Let  I7(fc)  be  the  L-point  DFT  of  u(n),  and  V{k)  the  /iT-point  DFT  of  v{n).  Then 

..  M-l 

t=0 

-  M-l 

=  ^  5]X(fc  +  i/f)E,„(fe-l-i/r) 

i=0 

„  X  M-l 

^^J2Xik  +  iK) 

i—0 

where  the  last  line  follows  because  we  have  Em{k)  =  Emik  -t-  iK).  But  the  last  line  is  merely  the  output  of 
the  system  shown  in  Fig.  10(b),  which  therefore  establishes  the  noble  identity.  Fig.  9(c)  shows  the  polyphase 
implementation  redrawn  using  this  result.  The  noble  identity  for  the  expander  is  similar. 
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A  caution  regarding  notation,  however,  is  in  order.  We  have  used  the  same  notation  Eni{k)  for  the 
mth  polyphase  component  in  Figs.  9(b)  and  9(c).  This  is  regarded  as  a  iiT-point  DFT  in  Fig.  9(c)  and  an 
L-point  DFT  (with  values  repeating  after  a  shorter  period  K)  in  Fig.  9(b).  The  (cyclic)  impulse  response 
of  the  filter  Em{k)  is  accordingly  as  shown  in  Fig.  10,  before  and  after  the  decimator.  If  this  distinction  is 
not  clear  from  the  context,  one  has  to  use  a  superscript  as  in  Em\k)  (before  decimator)  and  Em\k)  (after 
decimator). 

Now  consider  the  analysis/synthesis  system  of  Fig.  11(a).  With  the  filters  Hi{k)  represented  in  Type 
1  polyphase  form  and  the  filters  Fi{k)  in  Type  2  form,  we  have  the  equivalent  representation  of  Fig.  11(b) 
where  E(k)  and  R(fc)  are  the  polyphase  matrices  of  the  cyclic  filter  bank.  These  should  be  interpreted  as 
iif-point  DFTs,  e.g., 

K-l 

E(fc)=  (13) 

n=0 

The  filter  bank  has  the  perfect  reconstruction  property  x{n)  =  x(n)  if  and  only  if  the  equation  R(A:)E(A;)  =  I 
is  satisfied  for  all  the  K  values  of  k.  With  a  slightly  more  general  definition,  one  can  obtain  the  analog  of  Eq. 
(5.6.7)  in  [23].  The  alias-component  (AC)  matrix,  which  is  very  useful  in  the  noncyclic  case  [23],  can  also  be 
defined  for  the  cyclic  case,  and  the  alias-free  and  perfect  reconstruction  conditions  can  be  formulated  using 
this.  See  Appendix  B  which  also  shows  the  relation  between  the  cyclic- AC  matrix  H(A:)  and  the  polyphase 
matrix  E(A:),  which  can  be  used  to  express  orthonormality  directly  in  terms  of  the  filters  hi{n). 


3.  ALLPASS  AND  PARAUNITARY  PROPERTIES 


The  allpass  property  and  more  generally  the  paraunitary  property  play  a  crucial  role  in  digital  filtering,  and 
in  the  theory  and  implementation  of  multirate  filter  banks.  We  now  extend  these  ideas  to  the  cyclic  case. 


3.1.  Cyclic  Allpass  Filters 


A  cyclic(L)  allpass  system  is  one  for  which  H{k)  =  forO<fc<L>-l.A  simple  example  is  the  first 
order  system 


H{k)^ 


b*  +  Wl^ 


l  +  bWl 

We  can  always  rewrite  any  cyclic  allpass  system  H{k)  —  in  rational  form 


H{k)  = 


ELobnWt- 


C  =  1. 


(14) 


(15) 


For  example,  we  can  let  J  =  L  -  1  and  set  ^  The  coefficients  are  essentially  the 

inverse  DFT  coefficients  of  and  can  readily  be  identified.  More  interesting  is  the  open  problem  of 

obtaining  the  rational  form  with  smallest  order  J. 
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The  allpass  property  in  the  cyclic  case  can  be  expressed  entirely  in  terms  of  the  unitary  property  of  a 
circulant  matrix.  For  this  consider  the  L  x  L  circulant  matrix  H  formed  from  the  impulse  response 
demonstrated  below  for  X  =  4. 


H  = 


rft(O)  h{l)  h{2)  h{3) 
h{3)  /i(0)  h{l)  h(2) 
h{2)  h{3)  h{0)  h(l) 
./i(l)  h{2)  h(3)  /i(0) 


We  know  that  circulant  matrices  are  diagonalized  by  the  DFT  matrix  (which  is  unitary)  and  that  the 
eigenvalues  are  the  DFT  coefficients  H{k).  That  is, 

WAWt 


H 


M 


where  wtw  =  MI  and  A  is  a  diagonal  matrbc  with  elements  H{k).  Using  this  we  can  see  that  HtH  =  I  if 
and  only  if  |i/(A:)|  =  1.  Thus  the  cyclic  allpass  property  is  equivalent  to  the  unitariness  of  the  matrix  H. 

The  rational  form  (15)  has  nonzero  denominator  for  all  k  as  long  as  the  L-point  DFT  of  {6„}  is  nonzero 
for  all  k.  Equivalently,  the  polynomial  should  have  no  zeros  at  the  unit-circle  points  z  =  W^.  For 

the  allpass  case  we  can  assume  further  that  J2n  Kz~”'  has  no  zeros  anywhere  on  the  unit  circle  because  such 
a  zero  would  also  be  present  in  cancelled  anyway.  Thus,  any  cyclic  allpass  system 

cam  be  written  as  in  Eq.  (15)  where  bn  are  such  that  bnZ~'^  is  free  from  unit  circle  zeros. 


3.2.  Cyclic  Paraunitary  Systems  and  Orthonormai  Fiiter  Banks 

The  M  X  M  cyclic  transfer  matrix  E(fe)  is  said  to  be  cyclic-paraunitary  (or  cyclic-PU)  if  it  is  unitary  for  all 
k.  In  Sec.  4  and  5  we  consider  the  properties  of  cyclic-paraunitary  systems  in  greater  detail  and  show  that 
they  do  not  share  some  of  the  restrictions  of  noncyclic  PU  systems. 

We  define  the  M-band  cyclic  filter  bank  (Fig.  11(a))  to  be  orthonormai  if  the  polyphase  matrix  E(fc) 
is  unitary  for  all  k.  That  is,  E(A;)  is  cyclic-PU.  The  perfect  reconstruction  property  then  reduces  to  R{k)  = 
E^(fe),  which  can  be  rewritten  in  terms  of  the  cyclic(X)  impulse  responses  as 


Mn)  =  /i*(-n).  (16) 

The  DFTs  are  correspondingly  related  as  Fi{k)  =  H*{k).  Recall  here  that  the  arguments  of  hi{.),  /,(.), 
Hi{.)  and  Fi{.)  are  interpreted  modulo  L.  The  above  definition  of  orthonormality  is  consistent  with  the 
statement  that  the  filter  bank  expands  the  cyclic(L)  signal  x{n)  using  an  orthonormai  basis.  Thus,  assuming 
i(n)  =  x{n)  in  Fig.  11(a)  we  have 

M-lK-l 

=  I]  I]  -  ^M),  0  <  n  <  L  -  1.  (17) 

i=o  e=o 
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The  basis  functions  are  the  length-L  sequences 


»?i.£(n)=/i(n  -  £M),  0<n<L-l  (18) 

where  Thus  there  are  MK  =  L  basis  functions  7?i,«(n).  It  can  be  shown  that 

the  unitarity  of  R(fc)  is  equivalent  to  orthonormality  of  the  basis  r]i,e{n)  (Appendix  B).  This  orthonormality 
can  be  reexpressed  as 

E  /i(«)/m(n  -  Me)  =  S{i  -  m)5{e)  (19) 

n=0 

As  in  traditional  filter  banks,  orthonormality  of  the  cyclic  filter  bank  implies: 

1.  Unit-energy  property:  \Mn)f  =  1  ^r  all  i  (and  similarly  for  hi{n)). 

2.  Power  complementary  property:  Eilo  ^  l-^i(*)l^  =  M,  for  all  k  (and  similarly  for  Hi{k)). 


4.  CYCLIC  VERSUS  NONCYCLIC  SYSTEMS 

There  are  several  basic  differences  between  the  behaviors  of  cyclic  and  noncyclic  LTI  systems.  To  demon¬ 
strate,  consider  the  determinant  of  a  paraunitary  matrix.  For  the  noncyclic  case,  this  can  be  shown  to 
be  an  allpass  function  [23].  If  E(fc)  is  cyclic  paraunitary,  then  the  same  result  can  be  proved,  that  is, 
[det  E(fc)]  =  =  allpass.  However,  a  difference  in  behavior  arises  when  we  try  to  relate  the  degree  of 

determinant  to  the  degree  of  the  system.  The  degree  (or  McMillan  degree)  of  a  noncyclic  system  lEnoni^)  is 
defined  as  the  minimum  number  of  delay  elements  required  to  implement  it.  By  analogy,  in  the  cyclic 
case  suppose  we  define  the  degree  of  E(fc)  to  be  the  minimum  number  of  cyclic  delay  elements  Wj^  required 
to  implement  E(fe).t  For  noncyclic  FIR  paraunitary  systems  the  degree  of  [det  E„on(2)]  is  equal  to  the 
degree  of  Enon(-z)  (see  [23])  but  the  same  is  not  true  in  the  cyclic  case.  For  example,  consider  the  cyclic 
paraunitary  system 

E(fc)  = 

Here  [det  E(ifc)]  =  1.  Thus,  regardless  of  the  degree  of  E(fc),  the  determinant  has  degree  equal  to  zero! 

Another  difference  pertains  to  factorizability.  It  is  well-known  that  noncyclic  FIR  paraunitary  systems 
can  be  factored  [23]  in  terms  of  degree-one  FIR  building  blocks.  But  in  the  cyclic  case  such  factorization  is 
not  always  possible  as  explained  at  the  end  of  Sec.  5.4. 

4.1.  The  Noncyclic  Counterpart 

t  Notice  that  the  cyclic  nature  of  time  makes  this  rather  unnatural.  For  example,  the  system  H{k)  = 
requires  L - 1  cyclic  delays  but  can  be  rewritten  as  H{k)  =  Wi'^  requiring  only  one  cyclic  advance 

operator 


cos  6{k)  sin0(fc) 
—  sin  d{k)  cos  d{k) 


0<k<L-l 


(20) 
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In  the  cyclic(L)  case,  any  transfer  function  can  be  expressed  in  the  form  H{k)  =  Yln=o 
noncyclic  counterpart  of  this  is  defined  as 

Hnc{z)  =  h{n)z-^.  (21) 

n=0 

This  can  be  regarded  as  an  interpolated  version  in  the  frequency  domain,  with  H{k)  representing  the  samples 
of  Hnciz)  at  the  unit-circle  points  z  =  W['^  =  Similarly  the  noncyclic  counterpart  of  the  M  x  M 

cychc(L)  system  E(A:)  =  e(ra)W^"  is  defined  by  E„c(z)  =  5Zn=o 

The  interpolated  version  or  interpolant,  however,  is  not  unique.  For  example,  we  can  find  a  noncyclic 
interpolant  Gi„t(z)  =  ^ith  IV  >  £,  -  1  such  that  GintiW£'‘)  =  Hik).  As  another  example. 


consider  the  cyclic  system 


1  — 

=  — _ - _ 

“  1-aWt 

n=0  ^ 


If  we  replace  with  z~^  to  obtain  an  interpolant,  then  the  answer  obtained  from  the  left  hand  side 
is  different  fi-om  what  we  get  firom  the  right  hand  side.  Two  possible  interpolants  in  this  case  are  there¬ 
fore  y^Zn  a”z“”  and  .  It  is  clear  that  the  noncyclic  counterpart  is  only  one  of  the  many  possible 

interpolants^. 

If  E„c(z)  is  paraunitary,  it  readily  follows  that  E(fc)  is  cyclic-paraunitary  because  each  k  corresponds 
to  a  special  z  on  the  unit  circle.  But  the  converse  does  not  hold  as  we  shall  see.  Thus,  cyclic  paraunitariness 
is  less  of  a  constraint  on  the  coefficients  e(n)  than  traditional  paraunitariness.  To  demonstrate,  consider  the 
second  order  cyclic(4)  transfer  function 

G{k)  =  0.5  +  0.5(j  -  1)^4''  + 


Using  the  facts  that  =  1  and  =  —1,  it  is  readily  verified  that  |G(A:)p  =  1  for  all  k.  Thus  G{k)  is 
aUpass  in  the  cyclic(4)  sense.  However,  the  noncyclic  counterpart  Gnci^)  —  0.b+0.5{j  —  l)z~^  +0.5jz~^  which 
is  an  FIR  filter,  is  evidently  not  allpass.  If  we  now  construct  the  M  x  M  polyphase  matrix  E(A:)  =  G(A:)Im 
(for  arbitrary  M)  then  it  is  cyclic  paraunitary,  though  the  noncyclic  counterpart  E„c(z)  =  Gnc{z)l  is  not 
paraunitary. 

As  a  second  example,  consider  the  cyclic(3)  analysis  bank 

’Hoik)] 

H(fc)=  Hiik)  =ao  +  aiW3*=  +  a2W|'' 

,i?2(fc). 

t  The  subscripts  nc  will  be  for  “noncyclic  counterpart”.  We  will  use  the  subscript  int  to  refer  to  any 
interpolated  version,  and  the  subscript  non  to  indicate  any  noncyclic  system. 
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where  ai  are  the  column- vectors  given  by 


1 

'l' 

1 

[‘1 

1 

1 

0 

1 - 

1  1 
to  >-• 
1 _ 

’  ^/Io 

0 

It  can  be  verified  that  Ht(A:)H(A:)  =  1,  that  is,  the  three  transfer  functions  Ho{k),Hi{k)  and  H2{k)  satisfy 
the  cyclic  power  complementary  property: 

\Ho{k)f  +  \Hiik)\^  +  \H2ik)\^  =  l,  k  =  0,1,2. 

Consider  the  noncyclic  counterpart  Hnc(-2^)  =  ao  +  By  explicit  calculation  we  find  that 

H„c(2)Hnc(2)  7^  1  (e.g.,  the  coefficient  of  which  is  aja2,  is  nonzero).  Thus  the  noncyclic  counterpart  is 
not  power  complementary,  though  the  cyclic(3)  system  is. 

4.2.  Nyquist  Property,  Linear-Phase  and  CQF  Design 

The  next  example  brings  out  another  difference  between  cyclic  and  noncyclic  filters.  Consider  a  cyclic(6) 
transfer  function 

Hoik)  =  (l  +  -  Wt  +  Wt  +  Wi'‘)/VE, 

whose  impulse  response  is  shown  in  Fig.  12(a).  We  can  regard  this  as  FIR  in  the  sense  that  /i(n)  is  nonzero 
only  on  a  subset  of  points  inO<n<I/  —  1.  (In  the  cyclic  case  this  is  the  only  FIR-definition  that  makes 
sense).  The  symmetry  of  the  impulse  response  implies  the  linear-phase  property,  as  explicitly  seen  from 

Hoik)  =  -|-2cos(7rA;/3)  +  2cos(27rA:/3)^/-\/5 

Now  consider  the  magnitude-squared  function  G(A:)  =  |Ho(fc)|^,  whose  cyclic(6)  impulse  response  gin) 
is  the  cyclic  convolution  of  hoin)  with  hoi-n).  By  exphcit  computation  one  verifies  that  p(0)  =  1  and 
fl(2)  =  p(4)  =  0  (Fig.  12(b)).  That  is,  p(2n)  =  Sin).  So  the  filter  Gik)  has  the  half-band  property  (in  the 
cyclic(6)  sense),  and 

G(fc) -t- G(3  +  fc)  =  2,  0<fc<5.  (23) 

Equivalently  H ik)  is  power  S3Tnmetric  [23]  in  the  cyclic(6)  sense: 

\Hoik)f  +  |ifo(3  +  fc)l^  -  2,  0  <  fc  <  5.  (24) 

That  is,  we  have  found  a  cyclic(6)  filter  Hoik)  which  is  both  linear-phase  and  power-symmetric.  This  is  not 
possible  [23]  for  the  noncyclic  FIR  case. 
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Using  the  above  example  we  can  construct  a  two-channel  cyclic(6)  orthonormal  filter  bank  where  the 
filters  are  nontrivial  linear  phase  filters.  (Such  constructions  are  not  possible  in  the  noncyclic  FIR  case  [23].) 
For  this  choose  hQ{n)  as  above,  and  the  remaining  three  filters  as 

hi(n)  =  {-irhlil  ^  n),  /fc(n)  =  /ij(-n)  (25) 

for  fc  ==  0, 1,  with  all  arguments  interpreted  modulo-6.  That  is,  Hx{k)  =  —Wj^HQ{k  —  K)  where  K  =  L/2, 
and  Fi{k)  =  H*{k),  The  preceding  is  an  example  of  a  cyclic(6)  version  of  the  CQF  design  of  Smith  and 
Barnwell  [18]. 

More  on  the  CQF  Design 

Given  a  cyclic(L)  impulse  response  g{n)  with  the  property  that  G{k)  >  0,  we  can  construct  infinitely  many 
cyclic(L)  two  channel  orthonormal  filter  banks  by  choosing 

iyo(fc)  =  e^'^Wy^,  (26) 

and  the  remaining  three  filters  according  to  (25).  We  can  regard  i?o(fc)  ^  cyclic  spectral  factor  of  G(fc). 
Recall  that  for  noncyclic  FIR  filters,  the  usual  definition  of  spectral  factors  allows  only  finitely  many  phase 
responses,  and  the  spectral  factors  all  have  the  same  length  N  +  1.  In  Eq.  (26)  however,  the  phase  response 
</>(ifc)  of  the  spectral  factor  Ho{k)  is  arbitrary.  In  particular,  the  choice  <^(fc)  =  0  would  yield  linear  phase 
analysis  filters  Ho{k)  and  Hi{k)  with  good  frequency  responses  (if  G{k)  is  a  good  lowpass  filter).  However, 
even  if  G{k)  is  cyclic(L)  FIR  (in  the  sense  that  g{n)  =  0  for  iV  <  n  <  L  -  AT;  see  Fig.  13)  the  cyclic  spectral 
factor  Hoik)  may  not  be  FIR  (i.e.,  hoik)  could  be  nonzero  for  all  n).  An  open  question  here  is,  what  is  the 
most  general  form  of  the  phase  response  (^(fc)  which  ensures  that  Hoik)  is  FIR  (possibly  with  length  iV  -f  1)? 
A  related  question  is  for  a  given  G(A:),  what  is  the  minimum-length  linear-phase  spectral  factor? 

5.  THE  PARAUNITARY  INTERPOLATION  PROBLEM 

Given  an  arbitrary  cyclic(L)  paraunitary  system  H(fc),  that  is,  a  sequence  of  M  x  M  unitary  matrices 

H(0),H(1),...,H(L-1),  (27) 

can  we  always  find  a  paraunitary  interpolant?  That  is,  can  we  find  a  system  Hintie^^)  =  b(n)e"'^‘^^,  such 

that  its  samples  at  the  unit  circle  points  z  =  agree  with  H(fc)?  This  is  the  paraunitary  interpolation 
problem.  For  the  scalar  case  (M  =  1)  this  becomes  the  allpass  interpolation  problem.  For  the  matrix  case 
(arbitrary  M)  we  distinguish  between  the  FIR  interpolation  problem  where  Fiintiz)  = 
the  HR  case.  In  the  HR  case  we  again  distinguish  between  rational  interpolants  (where  each  element  in 
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Hi„t(z)  is  a  ratio  of  polynomials  in  z~^),  and  irrational  ones.  Notice  that  if  the  order  of  the  interpolant 
N  <L,  then  is  the  unique  noncyclic  counterpart  H„c(z)  defined  in  Sec.  4.1,  which  is  in  general  not 

paraunitary.  Things  become  more  interesting  when  N  >  L. 


5.1.  Scalar  Allpass  Interpolation 


Any  scalar  cyclic  allpass  filter  can  be  written  in  the  form  (15)  which  shows  that  there  exists  a  noncyclic 
allpass  interpolant  =  (cE;[=o  ^J-n^"”)/(En=o  where  |c|  =  1.  This  interpolant  is  HR,  and 

in  general  does  not  have  all  poles  inside  the  unit  circle.  For  example,  suppose  iV  =  1  and  consider  the  cyclic 
allpass  filter  H{k)  and  its  interpolant  Hint{.z)  given  below: 


b*  +  Wj^ 

1  +  bWj^ 


Hintiz) 


6*  +  z-i 

1  +  6Z-1 


(28) 


Let  |6|  >  1.  Then  the  cyclic  allpass  filter  is  still  well  defined  (because  1  +  bW^  0  for  any  k).  But  the 
interpolant  Hint{z)  represents  an  allpass  filter  with  a  pole  outside  the  unit  circle.  This  raises  the  following 
question:  suppose  we  allow  the  allpass  interpolant  to  be  of  higher  order  N  (even  possibly  N  >  L).  That  is, 
we  now  have 


(29) 


_ ’  Bint{z) 

Can  we  show  that  there  exists  a  polynomial  2n=o  ^nZ~'^  of  some  order  N  with  all  its  zeros  inside  the  unit 
circle,  such  that  H{k)  =  Hint{W£^)l  While  the  answer  to  this  is  not  clear,  the  following  example  gives 
evidence  to  the  contrary. 

Example  I.  Minimum  phase  interpolants.  Let  B{k)  =  1  +  bW^  with  6  >  1  and  let  Bint{z)  =  ^nZ~’^ 

be  an  interpolant  B{k)  =  Bint^Wh^)-  example  we  restrict  the  coefficients  bn  to  be  also  real.  We 

can  find  unlimited  number  of  choices  of  N  and  {6n}  satisfying  these  conditions.  But  for  even  L,  there  does 
not  exist  even  one  choice  such  that  Bi„t(z)  has  minimum  phase  (i.e.,  all  zeros  in  |z|  <  1).  This  is  proved  as 
follows:  since  the  coefficients  {6n}  are  real,  a  necessary  condition  for  Si„t(z)  to  have  minimum  phase  is  that 
■Bint(l)  and  Bi„f  (— 1)  have  the  same  sign  [2].  By  setting  fc  =  0  in  B{k)  =  BintiWl^)  we  had  B(0)  =  Bi„t(l). 
Similarly  B{L/2)  =  Bi„f(— 1).  Using  6  >  1  we  therefore  conclude 


B<nt(l)  =  B(0)  =  1  +  6  >  0,  Bi„t(-1)  =  B(L/2)  =  1  -  6  <  0 


Thus  Bi„t(l)  and  B<nt(— 1)  could  never  have  the  same  sign  no  matter  how  we  choose  N  and  {6n}-  Summar 
rizing,  B{k)  does  not  have  a  minimum  phase  interpolant  Bi„f(z).  However,  it  is  still  conceivable  that  the 
ratio  Bint  (z) /Bint  (z)  has  some  cancellations,  thereby  resulting  in  a  stable  allpass  interpolant.  Moreover,  we 
have  not  considered  the  possibihty  of  a  complex-coefficient  interpolant  (which  is  conceivable  even  for  real  b). 
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5.2.  FIR  Paraunitary  Interpolation 


Given  the  unitary  sequence  (27),  can  we  always  find  a  M  x  M  paraunitary  interpolant  restricted  to  be  FIR, 
i.e.,  of  the  form  answer  is  evidently  no,  because 

the  interpolant  has  to  be  an  FIR  allpass  function  (which  cannot  be  more  general  than  a  mere  delay).  So 
we  assume  M  1.  If  we  make  the  restriction  AT  <  L  the  coefficients  h(n)  are  simply  the  inverse  DFT 
coefficients  of  Hj„t(W£’*)  and  the  interpolant  is  the  noncyclic  counterpart  defined  in  Sec.  4.1.  This  may  not 
be  parauniteiry  as  seen  from  the  examples  of  Sec.  4.1.  More  generally,  suppose  we  allow  N  to  be  arbitrarily 
large  but  finite.  Does  this  allow  us  to  always  find  a  paraunitary  interpolant?  In  general  the  answer  is  still 
no,  as  we  shall  demonstrate.  For  this  we  first  review  a  well  know  result  for  noncyclic  2x2  causal  FIR 
paraunitary  matrices  [23]: 

Theorem  1.  Let  Hnoniz)  =  IDn=o  h(«)^“"  be  a  2  x  2  causal  FIR  system.  Then  it  is  paraunitary  if  and 
only  if  it  has  the  form 


Hnon(^)  — 


Ho{z) 

Hx{z) 


ei6z-^Hi{z) 

-ei^z-'^Hoiz) 


(30) 


where  Hq{z)Hq{z)  +  Hi{z)Hi{z)  =  1  (power  complementary  property),  6  is  real,  and  no  is  any  integer  large 


enough  to  ensure  causality. 


❖ 


A  proof  is  included  in  Appendix  C  for  completeness.  Consider  now  the  cyclic(L)  example 

H(fc)  = 


eia(fc)  0 


(31) 


0 

which  is  evidently  paraunitary  for  any  arbitrary  (real- valued)  a{k)  and  P{k).  Suppose  there  exists  a  noncyclic 
causal  FIR  paraunitary  interpolant  Hi„t(z).  If  this  interpolant  is  diagonal,  then  Ho(e-'")  =  and  a{k) 

has  to  be  of  the  form  —2‘nkN/L  plus  a  constant.  Thus  for  arbitrary  a{k)  and  P{k),  there  is  no  diagonal 
FIR  paraunitary  interpolant  Hi„t(z).  How  about  a  non  diagonal  interpolant?  Since  it  would  have  the  form 
(30),  a{k)  and  /?(fc)  cannot  have  arbitrary  combinations  of  values.  Only  those  combinations  which  satisfy 
the  property 

/?(fc)  =  — a(fc)  —  2nkno/L  +  5 


for  some  S  are  allowed.  This  shows  that  cyclic  paraunitary  systems  do  not  in  general  have  noncyclic  FIR 
paraunitary  interpolants.  Whenever  such  an  interpolant  does  exist,  it  can  be  factorized  into  the  form  [23] 

Hi„t(z)  =  Hi„t(l)  -  uiuj  -I-  z”^uiuj)  -  U2u|  -t-  z"^U2u|)  . . .  -  UjU^  -t-  Z"^UjU^ j  (32) 

where  u*  are  unit  norm  vectors.  By  replacing  with  we  obtain 

H{k)  =  H(0)  (l  -  uiut  W^uiul)  (l  -  U2uJ  +  Wj^U2nl)  . . .  (l  -  ujuJ  +  W^ujuJ)  (33) 
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which  is  a  factorization  of  the  cyclic  paraunitary  system  H(fe)  in  terms  of  the  building  blocks  -  UiuJ  + 
W^Ujut  Thus,  whenever  a  noncyclic  FIR  paraunitary  interpolant  exists,  the  cyclic  system  can  be  factored 
as  above.  Conversely,  if  the  cyclic  system  can  be  factored  as  in  (33),  then  replacing  with  z  ^  we  obtain 
a  FIR  paraunitary  interpolant.  Summarizing  we  have  proved: 

Theorem  2.  FIR  interpolants.  Let  H(k)  be  a  cyclic(L)  paraunitary  system,  (i.e.,  H(fc)  unitary  for 
0  <  fc  <  L  -  1).  Then  it  has  a  noncyclic  causal  FIR  paraunitary  interpolant  Hi„t(z)  if  and  only  if  H(fc)  can 
be  factorized  in  the  form  (33)  where  Uj  are  unit-norm  vectors.  0 


5.3.  HR  Paraunitary  Interpolation 

If  we  do  not  restrict  the  interpolant  to  be  FIR,  then  we  can  always  find  a  paraunitary  interpolant  for  the 
unitary  sequence  (27).  For  this  we  simply  define, 


Hi„t(e^")=H(fc), 


2Trk  .  27r(fc  +  1) 

- <  w  <  — - - 

L  -  L 


(34) 


for  0  <  jfc  <  X  -  1.  Then  the  sampled  values  Hj„t(VF£’*)  are  evidently  equal  to  H(A:).  The  interpolant 
Hi„t(e^")  is  a  piecewise  constant,  and  has  discontinuities  at  the  frequencies  27rfc/L.  It  is  therefore  not  a 
rational  function  in  e-’"  (i.e.,  the  elements  in  the  matrix  are  not  ratios  of  polynomials  in  z).  The  following 
theorem  asserts  that  we  can  always  construct  rational  solutions  in  the  HR  case.  Stability  of  the  interpolant, 
however,  is  not  asserted. 

Theorem  3.  HR  interpolants.  Let  H(0),  H(l), . . .  H(L  -  1)  be  a  sequence  of  M  x  M  unitary  matrices. 
Then  there  exists  a  causal  system  with  rational  transfer  matrix  Hi„t(z)  such  that  *')  =  H(A:).  <> 

Proof.  The  crucial  building  block  is  the  matrix 


Ufc(z)=I  -  uu^  +  z  ^Ffc(z)uu'l’ 


where  Ffc(z)  is  a  rational  allpass  filter  and  u  is  a  unit-norm  vector.  We  can  verify  that  Ufc(z)  is  paraunitary, 
that  is,  Ufc(z)U*;(z)  =  I.  Suppose  the  allpass  filter  Ffc(z)  is  chosen  such  that 

We  can  regard  Ffe(z)  as  a  rational  allpass  interpolant  (Sec.  5.1)  with  samples  at  z  =  as  specified 

above.  With  this  choice  of  Ffc(z),  the  matrix  Ufc(z),  sampled  at  z  =  yields 

TT  ^  for  all  m  A: 

torm  =  k 
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Now  any  M  x  M  unitary  matrix  can  be  expressed  as  a  product  of  M  -  1  matrices  of  the  form  I  -  2uut. 
More  precisely  [23]  each  matrbc  H(A:)  in  the  given  unitary  sequence  can  be  expressed  as 

H(fc)  =  (l  -  2ui,fcuJ  fc)  (l  -  2u2,kA,k)  •  •  •  (l  -  2uM-i,fcuL-i,fc)  A(fc) 

where  are  unit-norm  vectors  amd  A(fc)  is  diagonal  with  nth  diagonal  element  We  can  find  a 

rational  allpass  filter  Fn(z)  such  that 

Then  the  diagonal  matrfac  Dfc(2)  with  diagonal  elements  F„(z)  has  the  unit-circle  samples 

)-\A(A:)  m  =  k. 

By  multiplying  matrices  of  the  form  TJkiz)  and  Dfc(z),  we  can  define  a  noncyclic  paraunitary  system  Gfe(z) 
such  that 

Gk{W[^  )-\H(fc)  m  =  k 

The  product  =  Go(2r)Gi(z) . . .  Gl-i(z)  then  represents  a  rational  HR  paraunitary  interpolant  for 

the  given  matrix  sequenence  {H(A:)}.  V  V  V 

5.4.  Summary  of  Interpolants  and  Factorizability 

It  is  well  known  that  noncyclic  FIR  paraunitary  systems  can  be  factored  [23]  in  terms  of  the  building  blocks 
I  -  Ujut  +  z'^Ujut,  where  uj m  =  1.  But  in  the  cyclic  case,  factorization  in  terms  of  (I  -  UjuJ  +  W^Ujut )  is 
not  always  possible.  In  fact  Theorem  2  shows  that  such  factorization  is  possible  if  and  only  if  there  exists  an 
FIR  paraunitary  interpolant.  However,  the  fact  that  there  always  exists  a  rational  HR  interpolant  (Theorem 
3)  means,  in  particular,  that  we  can  obtain  a  factorization  of  the  cyclic  system  H(A:)  in  terms  of  slightly 
modified  building  blocks.  These  have  the  form  Ufc(z)  and  Dfc(z)  discussed  above,  with  z~^  replaced  by 
everywhere. 


6.  STATE  SPACE  DESCRIPTIONS  FOR  CYCLIC  LTI  SYSTEMS 


In  Sec.  2.1  we  considered  cyclic  diflFerence  equations,  and  recursive  structures  for  cyclic  transfer  functions. 
State  space  descriptions  allow  us  to  generalize  these  ideas.  From  the  direct-form  structure  of  Fig.  6  we  can 
identify  a  set  of  N  state  variables  Vi{n)  (outputs  of  the  unit  delay  elements  W^)  and  obtain  equations  of  the 


form 


v(n  -H  1)  =  Av(n)  +  Bx(n) 


(35) 


y(n)  =  Cv(n)  -1-  Dx(n), 
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where  v(n)  is  the  state  vector  given  by 


v(n)  =  [vi{n)  V2{n) 


VN{n)f 


(36) 


More  generally,  given  any  structure  for  a  cyclic  LTI  system,  we  can  identify  the  outputs  of  the  cyclic  delay 
elements  Wj^  as  the  state  variables,  and  obtain  the  above  equations.  Since  this  system  can  have  multiple 
inputs  and  outputs,  we  have  used  bold  letters  x(n)  and  y  (ti)  above.  Repeated  use  of  the  first  equation  in  (35) 
yields  v(L)  =  A^v(O)  +  a  linear  combination  of  samples  of  x(n).  Since  all  the  time-indices  are  interpreted 
modulo-L,  we  have  v(L)  =  v(0),  and 


(I  —  A^)v(O)  =  linear  combination  of  samples  of  x(n). 


(37) 


Thus  we  can  identify  the  initial  state  v(0),  provided  I  -  A^  is  nonsingular,  i.e.,  no  eigenvalue  of  A  has  the 
form  Wj^  for  any  integer  tti.  In  other  words,  the  eigenvalues  of  A  should  not  be  at  the  unit-circle  points 
indicated  in  Fig.  2.  This  nonsingularity  condition  can  be  understood  in  another  way.  If  we  evaluate  the 
frequency  response  H(A;)  explicitly,  we  would  have  the  form 


H(jk)  =D  +  c(w7*I-a) 


P(fc) 


(38) 


det  (W^'^I  -  A) 

The  eigenvalue  condition  on  A  implies  that  the  preceding  denominator  det  (W^^l  —  A)  is  nonzero  for  all 
integers  k.  As  long  as  the  eigenvalue  condition  is  satisfied,  H(fe)  is  defined  for  all  k,  and  we  can  uniquely 
identify  an  “initial  state”  v(0)  for  any  input  sequence  {x(n)}. 

Even  though  the  expression  for  H(fe)  resembles  the  noncyclic  case  Hnon(^)  =  D  +  C(2l  —  A)“^B,  the 
impulse  response  matrix  h(n)  takes  a  slightly  different  form  as  compared  to  the  noncyclic  case  [23]: 

,  -1 

R  n  =  n 

(39) 


h(n)  = 


fD  +  CA^-^(l-A^)  ^B  n  =  0 
=  S  .  /  -\  -1 


CA”-^  ®  1  <  n  <  I,  -  1 


Notice,  for  example,  that  h(0)  7^  D,  which  is  a  departure  from  the  noncyclic  case.  These  differences  arise 
because  the  initial  condition  v(0)  is  predetermined  as  explained  earlier,  and  cannot  be  set  to  zero  (as  we 
would  in  the  noncyclic  case).  One  can  verify  by  explicit  computation  that  the  expressions  for  h(n)  and  H(fc) 
are  indeed  related  as  H(fe)  =  h(n)W£*^.  As  a  final  remark  we  mention  that  though  the  computation 

of  v(0)  in  general  requires  an  initial  overhead,  such  a  computation  followed  by  the  recursive  computation  of 
y(n)  as  in  Eq.  (35)  is  often  more  economic  than  direct  or  FFT-based  circular  convolution  of  x(n)  and  h(n). 

Similarity  transformations.  If  we  define  a  new  state  space  description  by  using  the  familiar  similarity 
transformation  Ai  =  T~^AT,Bi  =  T“^B,Ci  =  CT,  we  can  verify  that  the  new  system  (Ai,Bi,Ci,D) 
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has  the  same  h(n).  Reason:  we  can  verify  by  substitution  that  CA"  ^  B  is  unchanged  by 

the  similarity  transform  for  any  n  >  1.  Thus  we  can  find  equivalent  cyclic  state  space  realizations  by  using 
similarity  transforms.  Note  that  even  though  D  does  not  represent  h(0),  it  is  still  unchanged  in  the  similarity 
transformation. 


6.1.  Reachability  and  Observability 

The  ideas  of  reachability  and  observability  [6], [9], [23]  can  be  extended  to  cyclic  LTI  systems  but  there  are  some 
differences  from  the  traditional  noncyclic  case.  For  example  we  will  see  that  reachability  and  observability 


together  do  not  imply  minimality.  The  cyclic  LTI  system  is  said  to  be  reachable  if  we  can  arrive  at  any 
chosen  final  value  v/  for  the  state  vector  v(n)  at  any  chosen  time  n  by  proper  choice  of  the  input  sequence 
x(.).  To  quantify  this  consider  the  state  recursion  v(n  + 1)  =  Av(n)  +Bx(n)  again.  If  we  apply  this  L  times 
and  use  the  periodicity  conditions  v(n  +  L)  =  v(n)  and  x(n  +  L)  =  x(n)  we  find 


(I  -  A^)v(n) 


[B  AB  ...  A^-^B] 


rx(n-l)l 

x(n-2) 


Here  we  have  used  the  notation  that  for  any  positive  integer  i, 


Lx(n  —  L)  J 


(40) 


7Za,bO)=IB  AB  ...  A‘-iB]  (41) 

Let  N  denote  the  state  dimension  (size  of  v(n))  and  r  the  number  of  inputs  (size  of  x(n)).  Then 
is  &  N  X  ir  matrix  with  rank  <  N.  The  matrix  72-a,b(L),  in  particular,  has  size  N  x  Lr.  Assume  (I  —  A'^) 
is  nonsingular  for  reasons  explained  earlier.  It  is  then  clear  that  we  can  attain  any  value  for  the  state  v(n) 
at  any  time  n  by  application  of  a  suitable  input  x(n  -  l),x(n  —  2), . . .  ,x(n  -  L)  if  and  only  if  the  matrix 
TIa,b{L‘)  has  rank  N.  This  gives  a  test  for  reachability.  Now  two  cases  should  be  distinguished: 


1.  Let  N  <L.  Then  the  rank  of  TIa,b{N)  will  be  equal  to  that  of  11a,b{^)  (Cayley-Hamilton  theorem). 
Then  the  reachability  test  reduces  to  the  conventional  one.  Moreover,  in  the  nonreachable  case  we  can 
perform  the  usual  reduction  and  reduce  the  size  N  of  the  state  vector. 

2.  Let  N>L.  This  is  possible  in  the  mimo  case  (e.g.,  if  H(A:)  =  W^Ir,  then  N  =  r  regardless  of  L).  For 
this  case  two  subcases  are  possible: 

2.1.  The  rank  of  'R.a,b{L)  is  already  N,  so  the  system  is  reachable. 

2.2.  The  rank  of  UaM^)  i®  smaller  than  that  of  TZaM^)-  If  *1^®  I®  ®^I®®  I®®® 

perform  the  usual  reduction  and  reduce  the  size  N  of  the  state  vector.  If  the  rank  of  TIa,b{^)  I® 
already  N,  we  cannot  do  this,  but  we  might  still  be  able  to  perform  a  reduction  of  the  cyclic  state 
space  equations  as  we  shall  demonstrate  below. 


18 


State-observability  in  a  cyclic  LTI  system  can  also  be  defined  similar  to  the  traditional  case,  but  with  some 
subtle  distinctions  between  the  cases  N  <  L  and  N  >  L.  First  assume  N  <  L.  The  output  equation 
y(n)  =  Cv(ra)  -b  Dx(n)  can  be  repeatedly  applied  to  yield 


y(n) 

■  c  • 

y(n  -b  1) 

CA 

.y(n-biV-  1). 

.CA^-'. 

v(n)  -b  func.  of  x(n),  x(n  + 1), .  •  •  x(n  -b  iV  -  1) 


(42) 


Sc.aW 


The  initial  state  v(ti)  can  be  uniquely  found  from  the  N  samples  of  the  input  and  output  in  this  equation, 
as  long  as  the  matrix  Sc,a{N),  which  has  N  columns,  has  rank  N.  If  N  >  L,  the  preceding  equation  is 
not  meaningful  because  y(i)  and  x(t)  repeat  with  period  L.  In  this  case,  however,  we  have  a  very  unusual 
situation.  If  the  input  and  output  are  known  for  all  L  values  of  time,  then  in  particular  x(i)  is  known  for 
all  i  and  we  can  identify  the  state  v(n)  for  all  n  using  the  state  recursion.  Thus  the  notion  of  observability 


becomes  trivial  for  N  >  L. 


Example  2 

Consider  the  cyclic  system 

H{k)  =  1  +  aWl  +  a?Wf  -b  . . .  -b  (43) 


for  which  a  direct-form  implementation  is  shown  in  Fig.  14(a).  IVith  state  variables  as  indicated,  the  state 
space  description  (A,  B,  C,  D)  can  readily  be  identified,  yielding 


■0 

1 

0 

...  O' 

■O' 

0 

0 

1 

...  0 

0 

0 

0 

0 

...  1 

,  B  = 

0 

,  C  =  [a^-1  a^-2 

..o] 

D  =  l 

(44) 

.0 

0 

0 

...  0. 

.1. 

Note  that  the  number  of  state  variables  N  =  L-l.  From  the  preceding  we  verify  that 


•0 

0  . 

.  0 

1 

O' 

ra^-1  a^-2  . 

a 

0 

0  . 

.  1 

0 

0 

Sc,a{N)  = 

0  0^-1  . 

.  a2 

T^a,b{L)  = 

0 

1  . 

.  0 

0 

0 

.  0  0  . 

.1 

0  . 

.  0 

0 

0. 

(45) 


Since  N  =  L-l,  TIa,b{L)  has  size  (L  -  1)  x  L  and  Sc,a{N)  has  size  (L  -  1)  x  (L  -  1).  Both  of  these 
matrices  have  rank  N  =  L-l  (assuming,  of  course,  a  ^  0),  showing  that  the  structure  is  both  reachable 
and  observable.  Notice,  however,  that  the  system  H{k)  can  be  rewritten  in  the  recursive  form 


(46) 
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using  the  fact  that  Wj;  =  1.  This  yields  the  simpler  recursive  implementation  requiring  only  one  cyclic  delay 
(Fig.  14(b)).  We  can  verify  that  the  state  space  description  of  the  simplified  structure  is 

A  =  a,  B  =  l,  C  =  a(l-a^),  D  =  l-a^  (47) 


In  this  case  the  number  of  state  variables  iV  =  1.  One  readily  verifies  that  11a,b0-)  =  1  and  <5c,4(l)  =  «(!  - 
a^).  So  TlA,BiL)  and  Sc,a{N)  have  rank  N,  and  the  structure  is  reachable  and  observable  (assuming  a  0 
and  ^  1).  Thus  the  two  structures  shown  in  Fig.  14  are  two  reachable  and  observable  implementations  of 
H{k)  with  different  state  dimensions!  The  first  one  requires  L-  \  cyclic  delays  (the  Wj^  elements)  whereas 
the  second  structure  requires  only  one  cyclic  delay.  Notice  also  that  the  quantity  D  is  different  for  the  two 
structures  unlike  noncyclic  systems.  This  is  consistent  with  the  fact  that  for  cyclic  systems  h{0)  ^  D  but  is 
given  by  the  more  elaborate  expression  (39). 


Example  3 

Consider  the  2x2  cyclic  system  shown  in  Fig.  15(a),  and  assume  L  =  3.  The  number  of  state  variables  is 


N  =  A.  The  state  space  description  has 


o 

o 

o 

-0  0- 

A  = 

0  0  10 
0  0  0  1 

,  B  = 

1  1 

0  0 

o 

o 

o 

o 

.1  1. 

C  = 


0  0  10' 
10  0  0 


(48) 


Then  explicit  computation  shows  that 


■0 

0 

1 

1 

0 

0- 

■0 

0 

1 

1 

0 

0 

1 

1- 

-0 

0 

1 

0- 

'R-aMI')  = 

1 

0 

1 

0 

0 

1 

0 

1 

1 

0 

1 

0 

,  7^^.B(iV)  = 

1 

0 

1 

0 

0 

1 

0 

1 

1 

0 

1 

0 

0 

0 

0 

0 

t  <5c,>i(2)  = 

1 

0 

0 

0 

0 

0 

0 

1 

.1 

1 

0 

0 

0 

0. 

.1 

1 

0 

0 

0 

0 

0 

0. 

.0 

1 

0 

0. 

Thus  TZ,a,b{L)  has  rank  Z  <  N  which  shows  that  the  cyclic  system  is  not  reachable.  However,  TIaM^) 


has  rank  4.  Since  Sc,a{^)  has  rank  4,  so  does  Sc,a{N).  So  we  cannot  perform  state-reduction  using  classical 


techniques.  In  this  example,  however,  it  is  possible  to  perform  state  reduction  of  the  cyclic  system  by  simple 
manipulations  of  the  structure,  and  by  using  the  fact  that  =  1.  For  this  we  notice  the  identity 


1 

1 


1 

1 


=  (l+a;) 


1 

1 


which  shows  that  the  transfer  matrix  of  Fig.  15(a)  is  eventually 
H(*)  = 


1  1' 

1  +  W|* 

0 

1  1 

_w^  +  wi'‘_ 

[1  1] 


(49) 


which  has  the  implementation  shown  in  Fig.  15(b)  requiring  only  two  cyclic  delays.  Thus  in  this  example, 
'R'A,b{N)  and  Sc,a{N)  have  rank  N  but  '1Za,b{L)  does  not,  and  we  have  been  able  to  reduce  the  state 
dimension. 
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In  Example  2  we  found  that  the  state  dimension  could  be  reduced  even  though  the  cyclic  system  is 
reachable  as  well  as  observable.  In  Example  3  we  found  that  'R.a,b(^)  Sc,a{N)  have  rank  N  and 
'R,a,b{^‘)  has  deficient  rank,  and  the  state  dimension  could  again  be  reduced.  The  question  now  is,  what  is  a 
necessary  and  sufficient  condition  for  the  minimality  of  state  dimension  in  cyclic  LTI  structures?  A  related 
question  is,  can  we  develop  a  theory  paralleling  the  Smith-McMillan  form  and  relate  the  minimum  state 
dimension  (McMillan  degree)  to  this  form?  These  appear  to  be  fundamental  questions  requiring  future  work. 


6.2.  Unitariness  of  Reaiization  Matrix 

Suppose  we  are  given  an  implementation  for  a  cyclic  transfer  matrix  E{k).  This  implementation  has  a  state 
space  description  of  the  form  (35).  The  realization  matrix  for  the  implementation  is  defined  as 


A  B 
C  D 


(50) 


The  following  result  connects  the  cyclic-paraunitary  property  to  unitariness  of  the  realization  matrix. 
Lemma  1.  If  the  realization  matrix  is  unitary,  then  the  cyclic  system  E(fc)  is  paraunitary. 

Proof.  Rewrite  the  state  equations  as 


’v(n-l- 1) 

A  B 

'v(n)‘ 

y(n) 

C  D 

x(n)_ 

Unitariness  of  the  realization  matrbc  implies  ||v(rH-l)|p  +  ||y(n)p  =  ||v(n)lp  +  ||x(n)|p  where  ||v||2  denotes 
v^v.  If  we  write  the  preceding  equation  for  0  <  n  <  L  —  1  and  add  them  up,  we  obtain 

L-l  L-l 

^yt(n)y(n)  =  53xt(n)x(n) 

n=0  n=0 

by  using  the  fact  that  v(n-f  L)  =  v(n).  With  X(fc)  =  and  Y(A:)  =  X)n=o 

obtain  (using  Parseval’s  relation)  Yt(fc)Y(A:)  =  xt(fc)X(A:),  that  is, 

L-l  L-l 

^xt(fc)Et(A:)E(fc)X(A:)  =  ^xt(fc)X(fe) 
fc=0  fc=o 

This  should  hold  for  all  sequences  {X(A:)},  which  implies  that  xt(fc)Et(A:)E(fc)X(A:)  =  xt(fe)X(fe)  for  any 
X(fc),  proving  Et(A:)E(A:)  =  I  indeed.  V  VV 

This  result  is  analogous  to  a  result  in  the  noncyclic  case  [23].  However,  unlike  in  the  noncyclic  case, 
we  do  not  have  the  converse  result.  That  is,  even  if  E(A;)  is  paraunitary,  there  may  not  exist  a  minimal 
nonrecursive  structure  (i.e.,  minimal  structure  with  all  eigenvalues  of  A  equal  to  zero),  with  unitary  system 
matrix.  When  such  a  structure  does  exist,  the  FIR  interpolant  Eint(^)  =  D  +  C(zl  —  A)“^B,  obtained  by 
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replacing  with  z~^  in  the  structure,  would  be  paraunitary  (because  a  result  like  Lemma  1  also  holds 
in  the  noncyclic  case  [23]).  Since  FIR  paraunitary  interpolants  do  not  always  exist  (Sec.  5.2),  the  point  is 
proved.  By  combining  this  observation  and  Theorem  2  we  obtain  the  following  result. 

Theorem  4,  Let  E(fc)  be  cyclic  paraunitary.  Then  the  following  statements  are  equivalent. 

1.  There  exists  a  causal  FIR  paraunitary  interpolant  Eint{z)- 

2.  E(fc)  can  be  factorized  into  unitary  building  blocks  as  in  (33). 

3.  There  exists  a  cyclic  implementation  of  E{k)  such  that  the  realization  matrix  is  unitary.  0 

7.  CONCLUDING  REMARKS 

The  main  purpose  of  this  paper  has  been  to  introduce  the  idea  of  cyclic  LTI  systems,  place  in  evidence  some 
interesting  theoretical  properties,  and  point  out  a  few  open  problems.  The  emphasis  has  primarily  been  on 
cyclic  versions  of  recursive  difference  equations,  allpass  filters,  paraunitary  matrices,  multirate  filter  banks, 
and  state  space  theory.  It  will  be  interesting  to  figure  out  how  to  exploit  the  extra  freedom  offered  by  the 
cyclic  system,  e.g.,  for  the  design  of  subband  coders.  Can  this  be  exploited  to  obtain  increased  coding  gain 
(or  compression),  or  to  reduce  the  complexity  of  implementations?  Evidently  more  work  is  necessary  in  order 
to  assess  the  practical  advantages.  We  saw  that  cyclic  LTI  systems  open  up  interesting  problems  in  the  more 
general  area  of  signal  and  system  theory.  Many  of  these  were  mentioned  throughout  the  paper.  One  can 
also  attempt  to  formulate  cyclic  versions  of  other  standard  problems  in  filter  bank  theory,  for  example  cosine 
modulated  filter  bank  design.  The  important  thing  again  would  be  to  establish  a  clear  advantage  for  such 
extensions. 

Appendix  A.  Proof  of  cyclic-decimation  formula 

By  definition  the  output  of  the  cyclic  decimator  is 

y{n)  =  x(Mn)  = 

m=0  m=0 

.  K-1  M-1  ^  K-1  .  M-1 

=  t  E  E  + ‘)  =  K  E  E 

fc3=0  i=0  k=0  t=0 

So  the  /T-point  DFT  of  y(n)  is  X(Ki  +  k)/M,  or  equivalently  X{k  -  iK)/M. 

Appendix  B.  Cyclic  alias-component  matrix 

Using  the  decimator  and  expander  formulas  (7)  and  (9)  we  can  express  the  L-point  DFT  of  the  reconstructed 
signal  x(n)  as 

M-1  M-1  - 

X{k)  =  V  X(Jfc  -  iK)  y\  ^Hm{k  -  iK)Fm{k),  0  <  A:  <  L  -  1, 

^  ^  ^  J\d 

1=0  m=0 
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where  K  =  L/M.  The  alias  components  are  X{k  -  iK),i  0,  and  the  desired  signal  term  is  X{k).  So  the 
perfect  reconstruction  condition  is 

Ml 
0 

0  . 

where  H(fc)  is  the  alias  component  (AC)  matrix  with  (i,m)  element  [H(fc)]im  =  Hmik  —  iK).  Now  the 
analysis  filters  and  polyphase  matrix  are  related  as 


H(fc) 


Fi(fc) 


[Hoik)  Hiik)  ...  HM-iik)]  =  [l  Wt  ...  wi^-^^'‘]E^ik) 


Using  the  fact  that  E(fc)  is  periodic(ii’),  we  can  show  that  H(fc)  =  WJ^A(fc)E^(A;)  where  Wm  is  the 
M  X  M  DFT  matrix,  and  A(A:)  is  diagonal  with  ith  element  Wj^.  Thus  orthonormality  (unitarity  of  E(A:)) 
is  equivalent  to  Ht (fe)H(A:)  =  MI.  In  the  two  channel  case  the  cyclic  AC-matrix  is 

Hoik)  Hi{k) 

Ho{k-K)  Hi{k-K)_ 

where  K  =  L/2.  If  Hoik)  satisfies  the  cyclic  power  symmetry  property  |//o(^)P  +  \Hoik  —  K)\^  =  2,  and  if 
Hiik)  =  ^Wj^H^ik  —  K)  we  can  verify  by  substitution  that  Ht(A:)H(fc)  =  21,  and  the  cyclic  filter  bank  is 
indeed  orthonormal.  Finally,  the  theory  of  alias-free  filter  banks  (i.e.,  the  pseudocirculant  conditions  [23]) 
can  be  extended  to  the  cyclic  case  in  a  straightforward  manner,  so  we  skip  the  details. 


Cyclic  orthonormal  basis.  For  orthonormal  cyclic  filter  banks  with  perfect  reconstruction  R(A:)  is  also 
unitary,  so  the  matrix  F(fc)  whose  elements  are  [F(fc)]rm  =  Fmik  —  rK)  also  satisfies  Ft(A;)F(A:)  =  MI, 
that  is,  -  rK)F^ik  -  rK)fM  =  -  m).  The  left  hand  side  is  the  AT-point  DFT  of  the 

M-fold  decimated  version  of  9i^)=^n=o  fiMfmi'^  ”  ^)*  orthonormality  of  R(fc)  is  equivalent  to 

Zln=o  /i(^)/m(^  “  ”  m)Si£)  which  proves  (19). 


Appendix  C.  Proof  of  Theorem  1 

The  FIR  matrix  H(2:)  is  of  the  form  H(z)  = 

^  ^  ^  ^  [ifi(z)  Gi(z)_ 

on  the  “only  if”  part.  The  paraunitary  property  H(z)H(z)  =  I  implies,  among  other  things,  that 


.  Since  the  “if”  part  is  obvious,  we  concentrate 


Hoiz)Hoiz)  +  Hiiz)Hriz)  =  l  (51) 

Hoiz)Goiz)  +  #i(z)Gi(z)  =  0  (52) 


Suppose  an  element  in  H(z)  is  identically  zero,  e.g.,  Hiiz)  =  0  for  all  z.  Then  Hoiz)Goiz)  =  0  for  all  z,  from 
(52).  Since  Hoiz)  is  nonzero  firom  (51),  we  conclude  Go(z)  =  0  for  all  z.  Thus  the  FIR  paraunitary  matrix 
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H(z)  becomes  diagonal,  proving  the  desired  form  (30).  So  we  assume  that  none  of  the  four  elements  in 
H(z)  is  identically  zero.  Now,  any  pair  of  power  complementary  transfer  functions  (i.e.,  a  pair  Ho{z),Hi{z) 
satisfying  (51))  should  be  free  from  nontrivial  common  factors  (i.e.,  factors  other  than  cz"”®),  because  the 
left  side  of  (51)  should  be  nonzero  for  all  z.  So  Hoiz)  and  ffi(z)  have  no  nontrivial  common  factors.  The 
same  holds  for  the  pair  Go(z)  and  Gi(z)  and  the  pair  Ho{z)  and  Hi{z).  So  from  Eq.  (52)  we  conclude 

Go(z)  =  coz-"'>#i(z),  Gi(z)  =  ciz-”^#o(2) 

which  upon  substitution  into  Eq.  (52)  yields  (coz~’^°  -\-ciz~”'^'jHo{z)Hi{z)  =  0  proving  cqz  +ciz  =  0. 
Thus  Co  =  —Cl,  and  no  =  ni.  Since  the  row  elements  Ho{z)  and  Go(z)  are  also  power  complementary,  i.e., 
Ho(z)Hoiz)  +  Go(z)Go(z)  =  1,  we  get  Hoiz)Hoiz)  +  \co\'^Hi{z)Hiiz)  =  1.  Comparing  this  with  (51)  we 
conclude  |coP  =  1,  i.e.,  co  =  e-^®.  This  establishes  the  form  (30).  VV  V 
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